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ABSTRACT 



Theoretical predictions of the added mass and damping coef- 
ficient curves for ships performing small vertical oscillations in calm 
water are determined. Calculations are made for the following ships: 

1. Series 60 Block Coefficient 0.70 

2. Series 60 Block Coefficient 0.60 

3. Golovato's Surface Ship (Weinblum) 

Model 

A method for the inversion of a power series type transform 
which conformally maps a unit circle into cylinders and ship-like 
forms is shown. Using this method, the transform coefficients cor- 
responding to the sections of the ships considered are obtained. From 
these transform coefficients, plots of the sectional added mass and 
damping coefficients against frequency of vibration are derived using 
Porter’s analytical solution to the problem. Comparisons are made 
between these curves and Grim’s predictions. 

Finally, by virtue of strip theory, the two-dimensional 
values for each cross section are integrated over the length of the 
ship to obtain the ship’s added mass and damping coefficients at each 
frequency considered. No correction is made for three-dimensional 
effects. The resultant curves for each ship are compared with the 
corresponding published experimental results. 
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I. INTRODUCTION 

The naval architect is confronted with the problem of predict- 
ing the forces that will act on a given ship moving in a given seaway in 
order to design for these forces. It has long been recognized that the 
solution to this complex problem can only be attained by breaking it down 
into a number of sub-problems, each one simple enough to hopefully 
permit an analytical solution. 

One such problem mentioned in the preceding paragraph is 
that of a rigid ship made to oscillate vertically by an externally-applied 
harmonic force in water that is initially still. The linearized force 
equation of motion in this case is of the form: 

m(l+kg)y + b^y + c^y = F^sin(a>t+ a) (1) 

This is a second order linear differential equation with constant coef- 
ficients whose solution is well known. However, k , the added mass 

s 

coefficient of the ship, and b^, the ship's damping coefficient, are usually 
not known. The problem is thus reduced to that of determining these coef- 
ficients for any particular ship at all frequencies of interest. The product 
mkg is called the added mass of the ship and accounts for the component 
of the hydrodynamic force acting in phase with the heave acceleration. 

The other component of the hydrodynamic force is the damping force 

which is the product b y and acts in phase with velocity. 

s 

Efforts to study the behavior of this complex hydrodynamic force 
led to a number of experiments. In particular, experimental results have 
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been made available for Series 60 Block 0. 70 by Gerritsma and Beukelman 
[1], [2], for Series 60 Block 0. 60 by Gerritsma [ 3] , and for a ship model 
with mathematically-defined lines* by Golovato [ 4] . Calculations made in 
this work are for these ships. 

In parallel with these experimental studies, considerable effort 
has been done to theoretically predict these quantities. In this field, the 
problem is further simplified from a three-dimensional ship form to that 
of an infinitely long cylinder with a ship-like cross section. With the 
acceptance of the validity of "strip theory", work along these lines for two- 
dimensional forms gained even greater importance. The cylinders of 
interest are those mapped by conformal transformation of a unit circle by 



z = 



; + 




n=0 



a r-(2n+l) 

^2n + 1 ^ 



( 2 ) 



The coefficients of this transform determine the shape of the section. 

Lewis [ 5] initiated the idea of fitting ship sections with forms from a 

two-parameter family of more or less ship-like forms which he generated 

by using this transform with chosen values of a^^ and a^. Landweber and 

Macagno [ 6] extended the generality of the sections generated by using a 

third coefficient, a^. for a three-parameter family of forms thereby 

b 

making possible a wider variation of forms. 



> 1 ' 



This model shall be referred to as the Weinblum Model in this paper. 
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The work most often referred to, at present, when theoretical 
predictions for the added mass and damping coefficients of two-dimensional 
forms are desired, is that due to Grim [ 7] . Grim generated a compre- 
hensive set of curves from which these quantities can be estimated for 
a given ship section using the beam-draft ratio and section area coef- 
ficient as the entering arguments. This is equivalent to approximating 
the section by a Lewis form, i. e. , by a shape defined by a^ and a^ only. 
Since a number of shapes can be defined by just one set of be am -draft 
ratio and area coefficient, it appears that two parameters would in 
general be insufficient to describe a particular ship section. This ambi- 
guity leads one to question the correctness of Grim's predictions. Now 
available is Porter's [8] analysis to the hydrodynamic problem which, 
within the limits of linearized potential theory, is exact. Since Porter's 
solution does not restrict the number of transform coefficients, we have 
in fact something that promises to give more accurate predictions. The 
author seeks an answer by using both methods. 

Porter's analysis assumes that these transform coefficients 
are known. Hence the need for a method to determine these coefficients 
for a given cylinder shape is called for. The first known method of in- 
verting this conformal transformation is that due to Plant [9] . This 
was further improved by Porter. While this method proved satisfactory 
in many cases, it was evident that a more accurate and less time- 
consuming method would be desirable. 

In aeronautics, the problem of conformally mapping an arbi- 

C I rc I €. 

trary airfoil shape into a has been solved by 
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Theodorsen [ 10] . His exact analysis requires the evaluation of a non- 
linear definite integral for numerical results. Naiman [11] reported a 
procedure for the numerical evaluation of this integral. The author 
devised a method of determining the transform coefficients for an arbi- 
trary ship section by an application of Theodorsen's method. The pro- 
cedure is outlined in this paper. 



M, ' K 
« ll<« > 10 

t 



» < !•«: tfy t . • ioa* * » fO- i 

1,1 . , •Ksmi*V» 

, ,, . r .i'.!iii » 1 I V fn*Maiil8V» UsI’i'^'Wn -jrff i*uib«oo JJ 

, ,vl .. VO. iM i jW«e*tr tui: to bOfUwo a bvalVab 

i I WI.II*»M« r nr.K iuw i-iHT to .dv- !» i» r*oHa!»« V.'toll 

.iw,9if L-^tdhuo a) anubao 



5 



II. PROCEDURE 



2. I Problem Statement 



Consider a rigid ship floating in smooth water. The depth 
of the water is infinite and its lateral dimensions are likewise infinite. 

We now impress upon the ship a vertical simple harmonic force so 
that the ship oscillates up and down in simple harmonic motion with a 
small amplitude. Assume that steady-state conditions have been at- 
tained so that the amplitude of the outgoing waves generated by the 
ship motion remains constant with time at any point on the water surface. 

Required for any given frequency are: 

a. The distribution of the added mass and 
damping coefficients along the length of 
the ship. 

b. The ship's added mass and damping 
coefficients. 

2. 2 Solution to the Problem 



Since the problem as formulated is not solved, we approxi- 
mate a solution by "strip theory". We divide the ship into a convenient 
number of stations with end stations at the forward and after perpen- 
diculars. To be definite, we will call the station at the forward perpen- 
dicular "station 1" and number succeeding stations consecutively up to 
the last station at the after perpendicular. Since in a later process we 
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will integrate over the length using an arbitrarily chosen method of 
integration, namely, Simpson's First Rule, the number of stations must 
be odd and the station spacings must be equal. 

We now proceed to solve for the sectional added mass and 
damping coefficients. The conditions stated in the problem statement in 
Section 2. 1 apply except that we replace the ship with an infinitely -long 
cylinder whose lower half cross-section is that of the ship's section up to 
its free- floating draft at a particular station. The axis of the cylinder is 
on the free surface of the water. It is clear that we need to consider as 
many cylinders as there are stations. Two-dimensional conditions and 
potential flow are implied. 

We define an added mass coefficient, as the ratio of 

the added mass of the cylinder to the mass of the fluid displaced by a 
circular cylinder of equal beam. For a cylinder with a half-beam, b, at 
the free surface, Porter showed that 



, , _ 4 4 N^A 

^ ” a" 



(3) 



and the corresponding damping coefficient, referred to the same circular 
cylinder, is 



c 



2it 






(4) 



For a given section, A, B, M^, and are functions of only the non- 
dimensional frequency, 6 = w b/g, in which w is the circular frequency 
of oscillation and g is the acceleration of gravity. A numerical pro- 
cedure to evaluate these functions is outlined in detail in Porter's 
paper [8]. Since the coefficients of the transform (2) determine the 
shape of the given section, we note that the added mass and damping 
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coefficients are functions of the frequency and transform coefficients 
only. It is clear that the number of these transform coefficients 
fixes the value of N in (Z). Thus we have the solution to the problem 
provided we know these "a's". 

The determination of these transform coefficients consti- 
tutes a separate problem. Grim infers fitting with two parameters. 
Landweber and Macagno propose three parameters. We introduce 
the use of the method described in Appendix A-1 which gives us a more 
accurate fit by a proper choice of N. The theoretical value of N for 
any given shape at present is still undetermined. We therefore 
arbitrarily set a reasonable criterion for our "best fit", and hence for 
our selection of N. 

Knowing the beam-draft ratio and the area coefficient of the 
cylinder, we can solve for unique values of a^ and a^ (see Appendix 
A-Z). By applying these values to equations (3) and (4), we obtain 
Grim's predictions. We note therefore that Grim's method is a special 
case of Porter's general solution. 

We now have two sets of values for the sectional added mass 
and damping coefficients at any frequency for each ship's station, 
namely, those derived from our "best fit" and those due to the method 
of Grim. We compare these two sets of results which are the predicted 
curves for the longitudinal distribution of the added mass and damping 
coefficients. 

On the basis of "strip theory", the ship's added mass, mk , 
and damping coefficient, b , can be obtained by integrating over the 
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ship's length the corresponding two-dimensional quantities. For each 
station at each frequency, we first multiply the value obtained 

by pS where p is the density of the water and S = b . Likewise, we 

multiply the c value by pSoo. We then integrate the resulting values 
over the length of the ship using Simpson's First Rule to obtain the 
ship's added mass and damping coefficient at each frequency. We 
compare these curves with the published experimental results. 

2.3 Numerical Calculations 



Numerical results were obtained for the following ships: 

1. Series 60 Block 0.70 

2. Series 60 Block 0. 60 

3. Weinblum Model 

For uniformity, twenty one stations and twenty one waterlines from 
the keel to the designer's waterline were used for each ship. The 
offsets used for the Series 60 ships are those punched in computer 
cards available in the Department of Naval Architecture and Marine 
Engineering Library. The offsets used for the Weinblum Model are 
those generated by the IBM 7094 computer at the Computation Center, 
M. I.T. using the mathematical definition for the model's lines. All 
offsets are normalized with respect to the maximum half-beam and 
given up to the third decimal place. 

The section-fitting method allows one to choose the number 
of transform coefficients to use. In these calculations, no more than 
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five coefficients were used solely for the reason of saving on computer 
time. In each case where the improvement in the fit was not substantial, 
a lesser number of coefficients was used. The method was programmed 
by the author thereby allowing calculations to be done by the IBM 7094 
digital computer. The transform coefficients selected for each station 
of the three ships are shown in Tables I, II, and III. 

In solving for Grim's a^^ and a^, the parameter 

. _ draft 

half-beam 

was used instead of the beam-draft ratio. This obviously did not 
change the end results. The other parameter used was the section area 
coefficient. 



section area 
Z X half-beam x draft 

For the Series 60 ships, the values of c were derived from the published 
ones. In the case of the Weinblum Model, the section area coefficients 
were obtained by analytical integration. A computer program of 
Professor Porter which takes X and cr as input data was used to cal- 
culate aj^ and a^. The values of X, o-, a^, and a^ are presented in 
Tables IV, V, and VI. 

Numerical calculations for the sectional added mass and 
damping coefficients were likewise done by the IBM 7094 digital com- 
puter using Professor Porter's program. 

The data for the models used by the experimenters were also 
used in calculating the ships' added mass and damping coefficients in 
order that the resultant curves may be compared with the published 
experimental results. 
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TABLE I 





Transform Coefficients for Series 60 


Block 0. 


70 Ship 


Station 


^1 




^5 




"9 


1 

2 


58924 


-. 03769 


.00647 


. 00548 


— 


3 


28865 


-. 05512 


.00595 


— 


— 


4 


10604 


-. 06690 


.01227 


. 00152 


.00771 


5 


. 00286 


-. 08005 


.01031 


-.00186 


.00698 


6 


. 06121 


-. 09 564 


.00662 


- . 00458 


.005Z6 


7 


. 08877 


-. 11214 


.00370 


— 


— 


8 


. 09 604 


-. 12545 


-.00249 


-. 00091 


.00358 


9 


. 10543 


-. 13671 


-.00928 


. 00533 


— 


10 


. 10 543 


-. 13671 


-.00928 


. 00533 


— 


11 


. 10543 


-. 13671 


-.00928 


. 00533 


— 


12 


. 10413 


-. 13545 


-.00808 


. 00402 


— 


13 


. 10047 


-. 12610 


-.00391 


— 


— 


14 


. 09419 


-. 10551 


. 00486 


-. 00525 


— 


15 


. 08197 


-. 08109 


.01288 


-.00260 


— 


16 


.05187 


-.04343 


. 02686 


— 


— 


17 


00480 


. 00432 


.04346 


— 


— 


18 


10876 


. 06381 


.06217 


.00847 


— 


19 


-.31290 


. 11763 


.09148 


. 02231 


.00311 


20 


70043 


. 12650 


.09670 


.05817 


.0Z740 
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TABLE II 

Transform Coefficients for Series 60 Block 0. 60 Ship 



Station 


^1 


^3 




^7 




1 

2 


77186 


02877 


. 006 56 


. 00522 


— 


3 


55189 


03135 


. 00364 


— 


— 


4 


34603 


02787 


. 000 69 


. 00619 


— 


5 


18343 


03758 


. 00410 


. 00206 


.00658 


6 


06696 


05339 


. 00617 


- . 00165 


.00744 


7 


. 01148 


06992 


.00582 


-. 00620 


.00469 


8 


.05821 


08784 


. 00412 


-. 00838 


.00381 


9 


. 08669 


10264 


.00265 


-. 00641 


— 


10 


. 09 729 


11366 


.00043 


-. 00475 


— 


11 


. 10190 


12402 


-. 00566 


— 


— 


12 


. 09 832 


11822 


-.00126 


- . 00488 


— 


13 


. 09 603 


10342 


. 00293 


-. 00601 


— 


14 


. 08899 


-.07995 


.00919 


- . 00679 


— 


15 


. 07168 


-. 05008 


. 0 19 44 


— 


— 


16 


. 03123 


-. 00606 


. 03972 


-. 00038 


.00687 


17 


02575 


. 04386 


. 05626 


. 00212 


.00430 


18 


13091 


. 09 183 


.07244 


.01083 


— 


19 


32985 


. 11966 


. 09513 


.03366 


— 


20 


64328 


. 09 866 


. 09761 


. 05619 


. 02954 


21 


05123 


. 16732 


. 00725 


.01821 


— 



I 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 
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TABLE III 

Transform Coefficients for Weinblum Model 



^1 


^3 




^7 


58026 


04867 


00944 


.00639 


34470 


-.07217 


-.01064 


.00440 


19658 


08429 


-. 01010 


-.00024 


09651 


09191 


-. 00877 


-.00340 


02668 


09738 


00712 


-.00564 


02263 


10174 


-. 00552 


-.00700 


05703 


10531 


-. 00421 


-.00760 


07986 


10808 


-.00330 


-.00771 


09 29 5 


10988 


-. 00280 


-.00760 


09724 


11051 


-. 00264 


-.00752 



. 00360 
. 00450 
. 00476 
. 00460 
. 00424 
. 00383 
. 00352 
. 00339 
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TABLE IV 

Grim Parameters for Series 60 Block 0. 70 



Station 


X 


or 


^1 


1 


— 


— 


— 


2 


3.941 


. 801 


-. 59129 


3 


1. 814 


. 827 


-. 28221 


4 


1. 198 


. 855 


-. 08598 


5 


.955 


. 888 


. 02167 


6 


. 848 


.922 


. 07488 


7 


. 810 


.956 


. 09368 


8 


. 800 


.976 


. 09752 


9 


. 800 


.985 


. 09686 


10 


. 800 


.986 


. 09678 


11 


. 800 


.986 


. 09678 


12 


. 800 


.985 


. 09 686 


13 


. 800 


.9 80 


. 09723 


14 


. 800 


.9 63 


. 09847 


15 


. 811 


.9 29 


. 09508 


16 


. 844 


.978 


. 07416 


17 


.920 


. 808 


. 04131 


18 


1. 080 


. 723 


-. 03982 


* 

19 


1. 465 


. 619 


-. 20811 


20 


2. 857 


.493 


-. 55146 


21 


— 


— 


— 



* 






See text paragraph 



for alternate values. 



Ship 



00659 
02443 
04381 
-.06583 
-.08731 
-. 10908 
-. 12231 
-. 12829 
-. 12896 
-. 1289 6 
-. 12829 
-. 12496 
-. 11374 
-. 09138 
-. 12412 
-.01441 
. 04002 
. 10280 



14533 
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TABLE V 



Grim 


Parameters 


for Series 


60 Block 0. 


60 Ship 


Station 


\ 


cr 


^1 


^3 


1 

2 


7.921 


. 822 


-. 76856 


-. 00934 


3 


3. 509 


. 823 


-. 54729 


-.01642 


4 


2. 046 


. 807 


-. 33923 


-. 01218 


5 


1. 426 


. 827 


-. 17107 


-. 02583 


6 


1. 114 


.857 


-.05155 


-. 04572 


7 


.951 


. 896 


. 02322 


-.07057 


8 


. 862 


.926 


. 06739 


-. 09020 


9 


. 817 


.953 


.08983 


-. 10719 


10 


. 800 


.967 


.09818 


-. 11634 


11 


. 800 


.977 


. 09746 


-. 12288 


12 


. 800 


.973 


.09775 


-. 12026 


13 


. 800 


.955 


. 09912 


-. 10792 


14 


. 805 


.922 


.09876 


-. 08670 


15 


. 820 


. 864 


. 09418 


-. 04961 


16 


. 853 


.781 


. 07961 


. 00267 


17 


.932 


. 693 


. 03703 


. 05863 


18 


1. 103 


.600 


-. 05505 


. 11930 


19 


1. 49 5 


.489 


-. 23533 


. 18555 


20 


2. 59 7 


. 346 


-. 54751 


. 23302 


21 


1. 076 


. 582 


-. 04138 


. 13109 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 
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TABLE VI 

Grim Parameters for Weinblum Model 



X 


O ' 


^1 


^3 


4 . 211 


.905 


-. 58706 


-. 04723 


2 . 222 


.917 


-. 35199 


-. 07202 


1 . 569 


.9 27 


-. 20218 


-. 08672 


1 . 250 


.937 


-. 10043 


-. 09613 


1 . 067 


.944 


-.02896 


-. 10241 


.952 


.951 


.02179 


-. 10670 


. 879 


.956 


. 05728 


-. 10960 


. 833 


.959 


. 08077 


-. 11149 


. 808 


.961 


. 09 420 


-. 11254 


. 800 


.962 


. 09 857 


-. 11289 
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III. RESULTS 



3. 1 General 



It has been stated in Chapter II that the offsets used for the 
Series 60 ships are those that have been previously punched in com- 
puter cards. These offsets are tabulated in Appendix B (Tables XXI 
and XXII). On the other hand, the section area coefficients used are 
those published in David Taylor Model Basin Report No. 1712. It is 
to be made clear that the numerical results presented in this paper are 
based on these data unless otherwise specified. This clarification is 
necessary because it was discovered late in the process of this investi- 
gation that there were differences in some of the offsets used and those 
in the above-mentioned publication. 

3. 2 Selected Fits 



The resulting fits due to the selected transform coefficients 
are shown for all stations on the body plan of each ship in Figures 1, 2, 
and 3. The solid lines are the actual ship sections and the corresponding 
fits are in broken lines. Since the Weinblum Model is symmetrical about 
the midship section, only half of the body plan is shown. 

3. 3 Grim Fits 



The fits due to Grim's a^ and a^ are shown (broken lines) for 




I 
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representative sections (solid lines) of Series 60 Block 0. 70 ship only. 
Figure 4 shows the two-parameter fit for Station 19 which is a vee- 
shaped section and Figure 5 for Station 11, the midship section. Since 
the offsets for Station 19 were seen to be different from the published 
ones, the area coefficient of this section was recalculated by integration 
(using Simpson's Rule) of the offsets used. Based on this section area 
coefficient (cr = . 506), a new set of a^^ and a^ (a^^ = -.22185; a^ = . 17562) 
was derived for this station to determine the Grim fit shown in Figure 5. 

3.4 Fits to a Bulbous Section 



A special type of form found on the forebody of many ships is 
bulbous. Such a shape is certainly of great interest in an investigation of 
this kind. Since none of the three ships considered has a section of this 
form, the section at the forward perpendicular of the Mariner was chosen 
for study. A five-coefficient fit (dotted line) and a Grim fit (broken line) 
are both shown with the actual section (solid line) in Figure 6. The 
corresponding transform coefficients and Grim parameters are listed in 
Tables VII and VIII, respectively. 

3. 5 Longitudinal Distribution of Added Mass and Damping 



The calculated sectional added mass and damping coefficients 
are presented in tabular form for each station of the three ships at 
different values of a non-dimensional frequency, 
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6 

r 

in which b is the maximum half-beam of the ship. The sectional added 
mass coefficients due to the selected fits are shown in Tables IX, X, 
and XI and those due to the Grim Fits in Tables XII, XIII, and XIV. 

The sectional damping coefficients derived from the selected fits are 
shown in Tables XV, XVI, and XVII and those due to the Grim Fits in 
Tables XVIII, XIX, and XX. Values of 1<2^4 ^ plotted against 

6^ for typical sections of Series 60 Block 0. 70 ship in order to have a 
comparison between the two methods. Stations 11 and 19 were again 
chosen as the representative sections. Figure 7 shows the sectional 
added mass coefficient curves and Figure 8 the sectional damping coef- 
ficient curves. Curves derived from the Grim fits are labelled "Grim" 
and those due to the selected fits are labelled "Porter". The values of 
and c for the Station 19 Grim curves are based on the new set of 
aj^ and a^ calculated for this station. The numerical computations for 
all added mass and damping constants used Porter’s program, simulat- 
ing Grim's results for two-parameter fits and used in the general way 
for forms using more coefficients. 

3.6 Ship's Added Mass and Damping 



2 . 
to b 



m 



The theoretical predictions for the ship's added mass and 



damping coefficients are plotted against frequency on the same graph 
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with the corresponding experimental curve for comparison. These 
curves are labelled in the same manner as that in the preceding section. 
The experimental curves were plotted by using points picked off the 
published ones. 

In order to conform with the values obtained by experiment, 

the results due to the selected fits and the Grim fits for the Series 60 

ships were made dimensional. The Block 0.70 ship's added mass, 

mk (in kg-sec^/m) , and damping coefficient, b Gn kg-sec/m), are 
s s 

plotted against the circular frequency, w(in radians/sec), and shown in 
Figures 9 and 12, respectively. The corresponding quantities for the 
Block 0. 60 ship are shown in the same manner in Figures 10 and 13. 
Since the experimental results for the Block 0. 60 ship are available only 
for several values of forward speed and none for pure heaving motion, 
the curves given for the lowest speed (Froude number, Fn = 0. 15) 
were used. 

In the case of the Weinblum Model, the predicted curves 

are presented in non-dimensional form to agree with the published 

results. Figure 11 shows the predicted and experimental curves for 

the ship's added mass coefficient, k , plotted as a function of the non- 

s 

dimensional frequency 




Likewise, the corresponding values of the dimensionless ship's 
damping coefficient, 
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t 




b 



s 




where A is the displacement and L the length of the ship, are 

I 

plotted against 6 in Figure 14. 



Designer's Waterline 



Zl 




Figure K Selected Fits to Sections of Series 60 Block 0.70 Ship 



Designer's Waterline 



2Z 




Figure 2, Selected Fits to Sections of Series 60 Block 0.60 Ship 
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Designer's Waterline 




(L 



Actual Section 



Selected Fit 



Figure 3. Selected Fits to Sections of Weinblum Model 




I 
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Actual Section 

Grim Fit 
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Figure 6. Fits to a Bulbous Section (Section 
1 of the Mariner) 
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TABLE VII 



Transform Coefficients for Station 1 of the Mariner 



^1 


^3 


^5 


^7 


^9 


85868 


11785 


1892 


. 00633 


00054 



TABLE VIII 



Grim Parameters for Station 1 of the Mariner 






cr 



a 



3 



133. 333 



7. 407 



85971 



12730 
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IV. DISCUSSION OF RESULTS 



4. 1 Fits to Sections 



It is evident from Figures I, 2, 3, and 6 that the section- 
fitting method introduced in this paper gives very close fits to all 
the forms considered. We note that the accuracy of the fit is greater 
for sections that are somewhat elliptical as in the case of the sections 
of the Weinblum Model and the sections in the middle body of the 
Series 60 ships. However, the loss of accuracy is small for fits to 
more general shapes such as those at the bow and stern of the 
Series 60 ships. The inaccuracies may be largely attributed to the 
nature of the numerical procedures employed and possibly due to the 
limited number of transform coefficients used. 

Figures 4, 5, and 6 show that the fits due to the method of 
Grim are not as close as those due to the new procedure used in this 
paper. We observe in Figure 4 that the slope at the waterline of the 
actual section is inclined at an angle of aboutforty-five degrees from 
the waterline while that of the Grim fit is vertical. Since it is known 
that Grim*s method is restricted to Lewis forms, this result is to be 
expected. On the other hand, the slope of the corresponding selected 
fit, as shown in Figure 1, is more or less the same as that of the 

* 

actual section near the vaterline and fits as well as at any other point. 

Comparing Figures 4 and 5, we note that a two-parameter 
fit to the midship section of Series 60 Block 0. 70 ship is relatively a 

better fit than a two-parameter fit to Station 19 of the same ship. 

------------ 

The slope at the vaterline of the selected fit is also infinite but 
for infinitesimal draft. 
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This is again to be anticipated since a closer approximation with a 
Lewis shape can be done to the midship section than to Station 19. 

The fits to the Mariner's bulbous bow section in Figure 6 
are quite interesting. It is to be noted that while the Grim fit is not 
as close a fit as that of the corresponding selected fit, it nevertheless 
is surprisingly a good one. This result is rather unexpected since 
such a form has heretofore been considered to be beyond the limits 
of even the three-parameter family of forms of Landweber and 
Macagno. It seems, therefore, that more work need be done in this 
field. Since any further investigation is beyond the scope of the 
present work, we leave this worthy endeavor to future research. 

4. 2 Sectional Added Mass and Damping 



It has been pointed out that Porter's solution to the problem 
was used to calculate the values for the sectional added mass and 
damping coefficients from the selected transform coefficients and 
Grim's parameters. We bring up the question as to how many 
transform coefficients or parameters are necessary to adequately 
describe the shape of the cylinder considered. It had been noted that 
two-parameters may give a reasonable approximation to the form 
provided the actual section is somewhat of a Lewis shape. However, 
as may be seen from Tables I, II. and III, at least three transform 
coefficients were selected for all the sections considered. This is 
due to the fact that relatively much better fits are obtained by using 
more than two "a's". In cases where two transform coefficients might 



I 
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have been considered to give reasonably good fits, it was observed 
that those "a*s" did not differ much from Grim's a^^ and a^* The 
Weinblum Model bore this conclusion out rather convincingly. 

As we might therefore have anticipated. Figure 7 and 8 
show that the values of the added mass and damping coefficients due 
to the two methods for the midship section (Station 11) of Series 60 
Block 0.70 ship are not very different. On the other hand, bearing 
in mind the fact that Grim's aj^ and a^ do not very well fit Station 19, 
one would expect a pronounced difference in the curves for this 
section. However, such an expectation is not fully realized as shown 
by the plots of these values. It seems that a more thorough investi- 
gation is necessary and we will not make any conclusions that might 
just be premature. 

4. 3 Ship's Added Mass and Damping 



Comparing the curves of the predicted values of the ship's 
added mass and damping due to the two methods, we note that the 
difference is small for every case. We also observe that these 
theoretical curves are in good agreement with the corresponding 
experimental results. This is especially so for the Series 60 Block 
0.70 ship as may be seen in Figures 9 and 12. Since the experiment 
was performed quite recently, it is reasonable to assume that the 
values obtained are more accurate than those for the other two ships 
due to improved techniques and better instrumentation. These en- 
couraging results indeed reaffirm the correctness and practicability 
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of Porter's solution. 

A comparative study of Grim's method and its results is 
not quite reassuring. We have noted from the Grim fits that two 
parameters may give a very different form from the actual section. 
Hence one has reason to feel unsafe when using the procedure. 
Nevertheless, the predicted curves are very similar to those due to 
the selected fits, at least for these three ships considered. 

On the other hand, we have seen how closely we were able 
to approximate the ship form by properly selecting the transform 
coefficients of the ship's sections. Hence we are more assured that 
we are making calculations for the correct ship by this more general 
application of Porter's solution. 
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V. CONCLUSIONS 



In this chapter, we summarize the conclusions stated or 
implied in the discussion of results (Chapter IV). 

We conclude that accurate fits to ship sections can be 
obtained by a procedure such as that introduced in this paper. We 
further conclude that in general, two parameters are not sufficient 
to adequately describe the shape of the section. 

The conclusion that one is more assured of correct results 
in calculating the added mass and damping coefficients of a heaving 
ship by a more general application of Porter's solution than that 
associated to Grim's procedure is certainly valid. 
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VI. RECOMMENDATIONS 

The results of this study show that more work needs to be 
done in certain specific aspects of this subject. 

The study on a bulbous section done in this work is just 
a brief beginning of what may be done in an investigation of the 
hydrodynamic properties of this important ship form. It is significant, 
however, in that a good fit to the section was obtained even with just 
two parameters. This seems to indicate that good fits could be ob- 
tained to a wide variety of bulbous forms. It would certainly be 
interesting to see how the added mass and damping coefficients of such 
forms would behave with respect to change of shape as well as to the 
frequency of oscillation. It is therefore recommended that further 
research along these lines be conducted. 

Likewise, a more detailed study on shapes with inclined 
slopes at the waterline is called for. To establish the effect of such 
slopes on the sectional added mass and damping coefficients would 
certainly be an important contribution. 

The effect of the goodness of fit on the estimation of the hydro- 
dynamic forces in sway motion and roll motion remains to be studied. 
Similar procedures as those used in this paper could be used in applying 



Porter's solution. 
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NOMENCLATURE 



transform coefficient; n = 0, 1 ... N 

half-beam of section 
damping coefficient of the ship 

dimensionless damping coefficient of the ship 

maximum half-beam of the ship 

sectional damping coefficient 
buoyant force coefficient of the ship 

acceleration of gravity 
sectional damping coefficient 

added mass coefficient of the ship 

mass of the ship 
time 

heave displacement 

heave velocity 

heave acceleration 
a complex variable 

amplitude of externally-applied vertical 
harmonic force 

Froude number 

length of the ship 

submerged area of a circular cylinder with 
half- be am, b 
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a 

6 

6 

] 

I 

6 

w 

P 

\ 

cr 

A 



phase angle 

a non-dimensional frequency 
a non-dimensional frequency 

a non-dimensional frequency 
circular frequency of oscillation 
mass density of water 
draft to half-beam ratio 
section area coefficient 
displacement of the ship 
a complex variable 
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APPENDIX A - DETAILS OF PROCEDURE 



A-1 A Method of Inverting a Conformal Transformation 



A-1. 1 The Problem 

Consider the two complex planes shown below 





where ^ = ipe and z = x + iy = ire Given M points on the 

ship section as shown in the z plane, determine the coefficients of 
the transform 



z = 




^2n+l 



^-(2n+l) 



which will conformally map the unit circle in the ^ plane into the 
ship section in the z plane, 



*Only the quadrant 0^0 < tt/ 2 is of interest. 
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A -1.2 Solution 



We separate the transform into its real and imaginary 
parts. For the j point of the M known offsets, we have 



N ^ 

X. = p sin 0. + X] (2n + l)0. 



J 



3 



n=0 



3 



(5) 



N ^ 

y. = p cos 0. + X! ~ 2n+l cos(^n+l)0j (6) 



3 



3 



n=0 



The " ratio, H, is a known constant and the following relation- 
ship 



H = 



N 

P + 51 (a 

n=0 



/^2n+l^ 
2n+l/P ^ 



N 

+ E (-1^' 



n=0 



(7) 



can be derived. Hence we have from (5), (6), and (7) a set of 2 M -I- 1 
simultaneous equations. The unknowns are p, a^^, and 0^ where 
i = l, 3, ... 2N-tl and j = 1 , 2, ... M. 

We first determine the angles, 0^. By an application of 
Theodorsen's method' , the angle in the ^ plane corresponding to any 
given point in the z plane can be calculated. This is done by using 

^ ' dj V 

an intermediate z plane (not shown) where z = ae^ h The 
transformation 

=1= The interested reader is referred to Reference 10 for the details 
of the method. 



60 



z = z + 



maps the ship section in the z plane into a curve in the z plane 

which may be expected not to differ greatly from a circle. and 

J 

"Yj can be solved from the following developed relationships: 



2 sin^y = P + < a ^ 



where 



2 sinh^ijj = -p+/p^ ^^a^ 



P ^ ^ 2a ^ ^ 2 ^ ^ 



Next, we wish to map the curve in the z plane into a 
circle in the ^ plane using the general transform 

00 






n=l 



z = ^e 



Theodorsen showed that 



”2- 



■^r 



( 0-0 ) 

J — d0. 



4j( 0) cot 
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If we know 4* as a function of 9, we can evaluate the integral by 
Naiman*s [ll] numerical procedure. Specifying the values of 4* at 2h 
equally spaced intervals in the range 0^0^ Zn, the method gives 

n 

e = / ^ -i- cot (4> , -4'i ) 

f — f n 2n '^-k ^k' 

k=l 

where the summation is for odd values of k only and 



However, we only know 4* as a function of y. Since y may be expected 
to differ but little from 6 , we take 4'( y) as a first approximation of 4'(9). 

A second approximation to the dependence between 4* and 9 is given by 
4i(0+e), This enables us to calculate a better estimate of e. Thus we see 
that an iteration process is necessary to determine both 4* and € as a 
function of 0 correctly. Knowing 6(9), we can get €(v) and therefore e . 



J 



This enables us to determine the angle. 



0 . = y. + € . 
J 1 3 



for each of the given offsets. 

Substitution of these now known angles, 9., into (5), (6), and 

3 

(7) leaves us N + 2 unknowns (p and d^) in the 2M + 1 simultaneous and now 
linear equations where 



d. 

1 



^ ^2n+l 

2n+l 

P 



n = 0, 



N. 
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We now solve this set of over -determined equations for these unknowns 
minimizing the sum of the squares of the residuals (or individual errors) 

t 

to get p = and d^^ = d^^ , i = 1, 3, ... 2N + 1. Since we have a choice 
over the value of N, we may choose N to be that value that gives the least 
total squared error. 

We are interested in the values of a. for p = 1. Hence it is 

i 

seen that p is in fact the scale factor and the desired normalized 
^o 

transform coefficients are 

I 

d. 

a. = -i i = 1, 3, . . . 2N + 1 
Po 

This numerical procedure has been programmed by the 
author in order that the calculations may be done by IBM 7094 digital 
computer at the Computation Center, M.I.T. This program has been 
submitted to the Department of Naval Architecture and Marine Engineer- 
ing Library of this school so that it may be available for general use. 
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A-2 Procedure to Derive Grim's and a^ 

It is known that Grim's method is restricted to Lewis 
cylinders only. For a Lewis shape, the following relationships hold 






1 - a. 



1 + a, 



+ a 



+ a 



3 



3 



TT 1 - aj - 3a^ 
4 (l+a^)*^ - a^"^ 



Knowing X. and cr, we can solve for unique values of a^ and a^. 
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APPENDIX B - ORIGINAL DATA 



The only purpose of this appendix is to include the offsets 
used in this work for the Series 60 ships. As pointed out in the text, 
some differences exist between these offsets and those published in 
DTMB Report No. 1712. 

The offsets for the Series 60 Block 0. 70 ship are tabulated 
in Table XXI. Likewise, the offsets for the Series 60 Block 0. 60 ship 
are listed in Table XXII. These values are normalized with respect to 
the maximum half-beam. The waterlines are numbered from 1 at the 

I 

designers waterline down to 21 at the keel. The values given for the 
waterlines are distances from the DWL and are in error. They should 
all be multiplied by 0. 8 to obtain the correct values. 

These offsets have been punched in computer cards and 
are available for general use at the Department of Naval Architecture 
and Marine Engineering Library, Massachusetts Institute of Technology. 
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